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Abstract 

The quantum charged rigid niembrane model, which is a higher derivative theory has been 
considered to explore its gauge symmetries using a recently developed first order formalism [1]. 
Hamiltonian analysis has been performed and the gauge symmetry of the model is identified as 
reparametrisation symmetry. First class constraints are shown to have a truncated Virasoro alge- 
braic structure. An exact correspondence between the higher derivative theory and the first order 
formalism has been shown from the point of view of equations of motion. 
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1 Introduction 



Higher derivative field theories are inseparable from modern day theoretical physics. Long ago physi- 
cists started considering Lagrangians with higher time derivatives O El H]. Initially they were in- 
troduced to avoid infinities appearing in the scattering amplitudes. But, due to their distinctive 
properties, HD(from now on the term "HD" will refer to "higher derivative") theories find its place in 
various context of physics e.g. electrodynamics [21 [3], non-local theories [5], relativistic particle model 
with curvature and torsion [BllTlll], string theory [9], supersymmetry [lOllllj noncommutativive theory 
[12] , dark energy physics [13l [HI [15] , cosmology [HI [T71 [18| , infiation theory [19] , brane world scenario 
[16] . supergravity [20l [21]. In quantum gravity, Stelle showed that adding higher derivative terms 
can ensure renormalizability |22] although it breaks unitarity. But a suitable choice of the coefficients 
of the higher derivative terms can lead to unitarity too [23]. People constructed f(R) gravity where 
higher curvature terms were added to Einstein-Hilbert action and opened a vast sector of research. 
For HD gravity, the list is huge. Interesting features appeared when higher derivative terms were 
added to study Higgs mechanism [24J. Also, people working in one of the most exciting fields of recent 
theoretical physics like AdS/CFT correspondence have considered HD theories [25 1 126 1 I^T j 128] which 
indicate the importance and relevance of considering HD theories. 

Existence of gauge symmetries in theories with higher derivatives can be an interesting domain 
to study. For theories with single derivatives only, there exists well established Dirac's method 
[291 130] l3Tl l32l 133] . But HD theories have some extra difficulties while performing canonical analysis 
and needed careful observation. Whereas, Ostrogradski's method for performing Hamiltonian analysis 
[34] specifically for HD theories can be useful , but with an extra burden of nontrivial definition of 
the momenta. For a long period the method was used in various sectors for higher derivatives theo- 
ries. This method was presumably first applied in the invariant regularization of gauge theories [35]. 
Other applications were done in various examples like equivalence theorems for spectrum changing 
transformations [36] . relativistic particle model[6l [7], ReggeTeitelboim type cosmology [18], geodetic 
brane cosmology [37], and recently for unambiguos quantization of nonabelian gauge theories [38]. 
Other than this, an inspired first order formalism exists in the literature where the HD fields are 
considered as independent fields and usual Hamiltonian analysis can be performed (along with a triv- 
ial definition of the momenta) [H [39]. For abstracting the gauge symmetries there exist a powerful 
method [401 HIl H2] but only for first order theories with no higher derivative terms. Recently, we 
provided a general method for abstracting gauge symmetries with higher derivative theories [H [39] 
which we referred to first order formalism. We obtained some peculiar result in gauge symmetries 
of HD theories. We took the relativistic particle model with curvature [6] and found that there are 
two independent PFCs (primary firstclass constraints) but with only one independent gauge symmetry, 
which is clearly contrary to the accepted result which states that the number of independent gauge 
symmetries is equal to number independent primary first-class constraints [331 141] . Surprisingly, there 
appears two gauge symmetries viz. diffeomorphism and W-symmetry when we considered the mass 
term to be zero [T]. These results inspired us to consider a thorough analysis of gauge symmetries of 
models with HD terms (especially with curvature terms). Such a model is Dirac's membrane model 
for the electron ^43] 144]. 

Theories with extrinsic curvatures are frequently studied especially in string theory. Although, the 
concept is not new but recent inclusion of these in some physically interesting models added an extra 
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urgency to revisit the symmetry features of this type of surfaces. Due to extrinsic curvature effects 
there appear geometrical frustration when nematic hquid crystals are constrained to a curved surface 
|45j . Whereas, graphene too can be considered as electronic membrane and its rippling generates spa- 
tially varying electrochemical potential that is proportional to the square of the local curvature |46j. 
These extrinsic curvature terms also appear in various brane world senariofTTl H8l HOj \50[ 151] , Re- 
cently, This concept of extrinsic curvature in membranes also have been incorporated for studying 
fluid dynamics [52]. Generally these surfaces come into the picture where we consider the evolution of 
a surface with a background metric. The lowest dimensional generalisation is a point particle evolving 
in spacetime with a background metric [53]. Applying this idea, in 1962 an extensible relativistic 
model of the electron was proposed by Dirac [43]. With spherical symmetry, the model was in stable 
equilibrium due to its surface tension. In this paper we shall investigate the gauge symmetries of an 
updated version of the Dirac's membrane model for the electron where extrinsic curvature terms of 
the world- volume were included as second order correction terms [44] . 

The paper is organised in the following manner. In section 2 we gave a general overview of 
higher derivative theories and their conversion to first order formalism. Construction of the gauge 
generator and the master equation for extracting independent gauge symmetries is introduced in this 
section. Section 3 comprises mainly of a very brief introduction to the model of quantum charged 
rigid membrane, since literature available for the model and its variants. Section 4 is purely new 
as our main work is concentrated here. In this section we derive the equation of motion from the 
variational principle and perform Hamiltonain analysis of the model. Section 5 is devoted to find out 
gauge symmetries. Interestingly, the first class constraints form truncated Virasoro algebra. In section 
6 we show the equivalence between the higher derivative and the first order formalism via matching 
the equation of motion. Finally, we conclude with section 7. 

2 Abstraction of gauge symmetries for higher derivative theories: 
a first order formaUsm 

A general form for HD Lagrangian is given b}0 



where x = Xn{n = 1, 2, • • • , i/) are the coordinates and ' means derivative with respect to time, i^-th 
order derivative of time is denoted by x'^'^\ 

In the first order formalism, we convert the Higher Derivative Lagrangian ([T]) into a first order La- 
grangian by defining the variables qn^a (a = !> 2, z/ — 1) as 




(1) 



Qn,a — 



qn,a-l, (a > 1) 



(2) 



Due to redefinition of the variables there emerges the following constraints 




(3) 



^for an extended version of this first order formalism please see [T] 
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which can be added to the HD Lagrangian via the Lagranges multiphers \n,p{P = — 1). 

Consequently, we can write down an auxiUary Lagrangian , 

^ i.Q.n,ai Qn,ai An,/?) = L {qn^i, qn,2 " " " , qn,v-l, qn,v-l) + ^ (^n,/? " (ln,l3-l) K,fi , (4) 

Considering the Lagrangian multiphers to be independent fields in addition to the fields (?n,o, we define 
momenta as 

dV dV 
Pn,a = jr: — , vr„,^ = — . (5) 

Having found out the primary constraints of the theory, we can write down the total Hamiltonian as 

Ht = He + Un,j3TTn,P + ^^n,/3^n,/3, (6) 

where Un^j^^Vn^p are Lagrange multipliers and iTn,i3i^n,i3 are primary constraints. So we can proceed 
to have all the secondary constraints by demanding time variation of the constraints as zero. After 
we have extracted all the constraints, we can move to distinguish the first class and second class 
constraints. Now, according to Dirac, the first class constraints generate gauge transformation. The 
second class constraints can be removed by introduction of Dirac brackets. Therefore, our theory is a 
first order theory with only first class constraints. To find out the gauge symmetries of the model we 
define the gauge generator as 

G = Y.ea^a. (7) 

a 

Here {^a} is the whole set of primary constraints. All the gauge parameters may not be 
independent. To identify all the independent gauge transformation we refer to the method developed 
in [m 132] and write down the master equation relating the Lagrange multipliers Aa^ and the gauge 
parameters 

SAa, =^-^a (Kai + Afc, Cb.aa, ) (8) 
0=^-6, {Vaa, + Afei Cb,aa, ) (9) 

Here the indices ai,bi... refer to the primary first class constraints while the indices 02, 62... corre- 
spond to the secondary first class constraints. The coefficients V^^ and C^^^ are the structure functions 
of the involutive algebra, defined as H 

{H^an,'^a}D = Vab^b 
{'^a,'^b}D = Cabc'^c (10) 

Due to the HD nature, a relation between the gauge transformations of the fields can be written as 

Sqn,a - ^Sqn,a-l = 0, (q > 1) (11) 

which may impose some extra condition on the gauge parameters. 



^from now on we have to use only Dirac brackets since we removed all second class constraints. Poissson brackets are 
denoted by { , } , whereas, { ,}d refers to Dirac brackets 
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Abstracting all independent gauge transformations, we can write gauge transformation of the basic 
fields as 

SeaQn,a = {qn,a,G}D (12) 

This completes our analysis of finding the gauge transformation for HD theories. 



3 Quantum charged rigid membrane 

In a background Minkowski spacetime r/^'', H consider the evolving surface S. The surface is described 
by the local coordinate of the background spacetime. The embedding function X^{^°') = is a 
function of the local coordinates of the world volume m, swept out by the surface. We consider the 
following effective action underlying the dynamics of the surface E|44]: 



S[X'^]= [ <fi{-aK + Pfe\A^), (13) 

J m 

where K = g°'^Kah being the extrinsic curvature |l| and a, /3 are constant related to the rigidity pa- 
rameter and form factor respectively. On the other hand, j"" which minimally couples the charged 
surface and the electromagnetic field [55], is a constant electric current density distributed over 
the world volume and is locally conserved on m with daj"" = 0. Variation of the action with respect 
to the embedding function X'^{^°-) leads to the equation of motion 

an = -^fn'^e^F^,. (14) 

The above equation (|14|) can be thought as a Lorentz force equation with TZ being the Gaussian 
curvature and F^i, = 2d[^A^] the electromagnetic field tensor. Under suitable choice of the embedding 
functions {X^^{T,9,{p) = (t(r), r(r), 6*, 99)) equation boils down to[ 



5 = 47r / dTL{r,r,r,i,i) (15) 



where the Lagrangian L, which is HD in nature is given by, 

L = -^ -{ft-rt) -2art- (16) 

So, Lagrangian (|16p will be our sole interest which is reparametrisation invariant under the parameter 
r. Promptly, we can write down the equation of motion for the HD Lagrangian: 

2 -2 / ofl2 ^2\2„2\ 



Vith = 0, 1, 2, 3 and a, 6 = 0, 1, 2 



gab is the worldvolume metric and e^a = ^^,a a-re tangent vectors to the worldvolume 
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4 Hamiltonian analysis 



Before we start the Hamiltonian analysis we need to convert the HD Lagrangian (jl6p to a first order 
lagrangian, named as the auxiliary lagrangian, by introduction of the new fields 

r = R 

i = T (18) 
So, we write down the auxiliary Lagrangian as H 

L' = -^{RT - RT) - 2arT - + Ai(i? - r) + \2{T - t) (19) 

Inclusion of new fields impose constraints 

R-r^Q, T-i^O (20) 

which are taken care of via the multipliers Ai and A2. Variation of L' with respect to r, R, t, T, Ai and 
A2 give rise to the following equation of motions: 

-f^{RT-RT)-2aT + ^ + Xi = (21) 

A2 = (23) 



^^nRT-RT)--^-j^Rj-^R-2ar-^ + X, = (24) 

R-r = (25) 
T-i = (26) 



25|) . (j26|) are obvious since they correspond to (j20[) . 



Before proceeding for Hamiltonian formulation, we identify the new phase space which is con- 
stituted of the variables are (r, H^), (t, Hf), (i?, H^j), (T, Hy), (Ai, H^J, (A2, HaJ. Here H^jm = 
are the momenta corresponding to which generically stands for the variables r, R,t,T, Xi, X2- We 
immediately obtain the primary constraints as listed bellow 



$1 


= Ur + Xi^O 






= Ht + A2 ~ 






2 




$3 




i 




ar 






:i 


$5 


= « 




$6 







^ consider A'^'^ = — R^, for convenience 
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The poisson brackets between the field variables are defined as: 



K.x'^} = {n,.,n,.} = o (28) 

With the aid of (j28p the non zero Poisson brackets between the primary constraints can be written 
down 

{^1,^3} = -j^T 
, 2Q;r „ 

{^1,^4} = J^R 

{$1,^5} = 1 

{$2, $6} = 1 (29) 
We can take the following combination of the constraints 

= i2$3 + r$4«0 (30) 
= $4-^<I>5-0 (31) 

so that the new set of primary constraints are $1, $2, *^*3, *1*4, *1*5, *&6- The complete algebra of primary 
constraints is now given by (only the nonzero brackets are listed), 

{$1,$5} = {^2,^6} = 1 (32) 

We can write Canonical Hamiltonian via Legendre transformation as 

H,an = 2arT + ^ X^R- MT. (33) 

r 

The total Hamiltonian is 

Ht = Hcan + Ai^i + A2«'2 + As^s + A4^4 + ^5^5 + Ae^G (34) 

Here Ai, A2, A3, A4, A5, Ag are the Lagrange multipliers which are arbitrary at this stage. Only those 
multipliers which are attached to the primary second-class constraints will be determined, others 
corresponding to primary first class constraints will remain undetermined (although they can be 
determined too via equation of motion). At this level, loosely speaking $3 and are first class 
constraints (this classification may be changed after we get the full list of constraints). These two 
may provide us two new secondary constraints and the list can still keep increasing until we get all 
the constraints. Now, we move towards extracting all constraints of this system. This can be done by 
demanding that Poisson brackets of the constraints with the total Hamiltonian (time evolution) of the 
constraints is zero. Preserving ^1,^2,^5,^6 in time solves the following multipliers respectively 

A, = 2aT-^ 

Ae = 
Ai = R 

A2 = T. (35) 
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Whereas, time conservation of the primary constraints and $4 leads to the secondary constraints 
^1 and ^'2 respectively given by 

*i = -2arT-— — + XiR + X2T^0 
r 

*2 = -2ar - ^ + A2 - --^i?^ ~ 36 

Before proceeding further we list below all the nonzero Poisson brackets of the secondary constraints 
^1, ^2 with other constraints: 



.2 



{$1,^1} = 2aT - 
{cl>5,*i} = -R 
{$6,^-1} = -T 

{^4,^2} = -^Ti?^ 

{$6,^2} = -1 (37) 

Now, time preservation of the secondary constraint ^1 gives identically = 0. And requirement of 
^2 = solves the Lagrange multiplier A4 = -^i?, with A = ^ - ^ and 5 = 



Prom the constraint algebra (j32p and (j37p one can clearly assert that there is only one first class 
constraint ^'■^ with seven other second class constraints $1, $2, *1'45 'I's, ^6) ^i) ^2- One point worth 
noting since there are odd number of second class constraints, it indicate there might be some other 
first class constraint to make the pair of second class constraints even. Judiciously, we can choose 
a combination = xj/]^ - Ai$i - A2$2 - A4<i>4 - A5<I>5 - Ae^e so that the pair (<^3,^i) becomes 
first-class. This completes our constraint classification. 

Having completed the constraint classification, its time to get rid of the unphysical sector (Ai, 
and (A2, 11^2) by imposing the primary second class constraints 4>i, $2, $5, ^"6 strongly zero. This can 
be done by replacing all Poisson brackets by Dirac brackets for rest of the calculations. Surprisingly, 
Dirac brackets between the basic fields remain same as their corresponding Poisson brackets. So, now 
our phase space is spanned by {r, 11^, t, IIj, i?, II/j, T, 11^}. For convenience of future calculations we 
rename the constraints as 

Fi = % = R<^3 + r^>4 w 
F2 = ^i- A4<I>4 w 

51 = $4«0 

52 = ^-2 = -Ht - 2ar - — 

r 



(38) 
(39) 
(40) 



0. (41) 



Here, Fi,F2 is the first class pair with Fi as primary first class constraint. So far we observed that 
in this theory, there is only one primary first class constraint with one undetermined multiplier which 
clearly indicate existence of gauge symmetry(s) in the system. In the next section we will extract the 
gauge symmetries of this quantum charged rigid membrane. 
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5 Gauge symmetry and Virasoro algebra 



To study gauge symmetry we need to remove all the second class constraint from the system by setting 
them strongly zero and performing Dirac bracket defined by 

{f,9}D = {f,9}- {f^Si}A-^'{Sj,g} (42) 
where f and g corresponds to the phase space variables or their functions. To compute A^-^ for the set 



of of second class constraints, we have {81,82} = — "'J^-i ■ So, we can compute the Dirac Brackets 
between the basic fields. The nonzero DBs are: 

{r,Ur}j, = 1 

{Ur,t}^ - 
{Ur,t}j^ 



2TR 

A 



B 

„ , 2arR Ar 



{t,T}n 

{t,'nT}D 



1 

r 

~2R 
{R,IIr]d = 1 

{^R,T}n = -| 

2arT A ar'^ {T'^ + R^) 

|n„nR|^ = ^^ + ^ 

2 

{TIr,TIt}j, = ^ (43) 

The generator of the gauge transformation is given by a linear combination of all first class con- 
straints, 

G = eiFi + £2^2 (44) 

where ei and €2 are gauge parameters. We need to find out whether these gauge parameters are 

independent or not. 

The Dirac brackets between the first class constraints are given by 

{Fi,Fj}D = -eijF2 ; i,j = 1,2 (45) 

Using a suggestive notation we rename the constraints Fi and F2 as 

Lo = Fi (46) 
Li = F2 (47) 
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We can easily identify a sort of truncated Virasoro algebra of the form 



{Lm, Ln}D = {m- n)Lm+n (48) 



with m = 0, n = 1 as proposed in |54j for HD cases. 

Now, using equations (flOl 1451) we compute the structure constraints as C122 = — 1 = — C212 and 
V12 = 1 (other structure constraints are zero). Exploiting the master equations ([9]) we find the the 
following relation between the gauge parameters 

ei = -Agea - h (49) 

and It is clear that we have only one independent gauge symmetry in this system which is supported 
by the fact that there is only one undetermined multiplier. We consider £2 to be independent and 
compute the gauge transformation of the fields 

8r = -e2R (50) 
5t = -e2T (51) 
5R = eiR (52) 

5T = eiT + e2^R (53) 

We can identify this gauge symmetry as reparametrisation symmetry in the following manner. 
Consider an infinitesimal transformation of r and t on the worldvolume as r — t- t + o". For some 
infinitesimal a, we can write 

6r = —ar 

5t = -at (54) 



Clearly, a comparison between (150ll5ip and both equations of (I54p shows that the reparametrisation 
parameter is given hy a = €2- Using ()54p we compute of Gauge variation of the Lagrangian (jl6p which 
simplifies to 

6L = ^{aL) (55) 
ar 

and ensure the invariance of the action under ([5 



6 Consistency check 

It would be worth to find out the Hamiltonian equations of motion which are given by 

f = R (56) 

i = T (57) 

R = A3R (58) 

t = -iii + |r (59) 
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Equations (j56p and (j57p are obvious as they arise as constraints at the Lagrangian level and agrees 
with (p5]) and (j26|) . Taking time derivative of ([50|) and (fST]) we get 




- esT 



e2R-e2R 



(60) 



(61) 



Using equation (^9]) alongwith ([58l [59]) the above equations ([60l [6T]) simplify to 




6R 



(62) 



(63) 



which is a direct verification for (jlip . Whereas, (|59p along with the trivial equation of motions (|56p 
and ()57p can be cast into the form so that it verify ()17p . This indeed is an important outcome of this 
analysis which agrees the validity of this first order formalism via matching the equation of motion at 
higher derivative and first order level. 

Taking gauge variation of the equation of ([58|) and using ([52]) we get 



which in turn verifies the first master equation ([8]). 

7 Discussion 

Studies in higher derivative field theories have been an intense field of research [21 [3l [3 [8]. Symmetry 
studies has always been interesting for theoreticians. We already have shown some result concerning 
inequality in number of independent first class constraints and number of independent gauge sym- 
metries for a relativistic particle model with curvature[lj. This mismatch inspired us a further study 
of some physically interesting model. Dirac's relativistic membrane model for the electron can be a 
candidate with future prospect in brane inspired cosmology |47] . 

In this paper we presented a fresh Hamiltonian analysis purely in a first order formalism where 
higher time derivatives are considered to be independent fields and the corresponding momenta are de- 
fined in the usual way. Gauge symmetries were analysed with a novel way by constructing the gauge 
generator and extracting the independent gauge parameter. Number of independent primary first 
class constraint exactly is in accord with number of independent gauge symmetries leading to no mis- 
match. Also the constraint structure is shown to obey truncated Virasoro algebra. Reparametrization 
parameters have been identified through a suitable transformation of the fields. 

The model continues to be in the highlight of recent interests like branes, cosmology and dark 
energy [471 148| [52| [55| I56j. Consideration of other variants of the model with more symmetries can be 
of utmost interest as future projects . 



5 A3 = ei 



(64) 
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